Hard pomeron enhancement of ultrahigh-energy neutrino-nucleon 

cross-sections 



A. Z. Gazizov and S. I. YanushF] and S. I. YanushQ 

B. I. Stepanov Institute of Physics of the National Academy of Sciences of Belarus, 
F. Skariny Ave. 68, 220072 Minsk, Belarus 

Abstract 

An unknown small- x behavior of nucleon structure functions gives appreciable uncertainties to 
high-energy neutrino-nucleon cross-sections. We construct structure functions using at small x 
Regge inspired description by A. Donnachie and P. V. Landshoff with soft and hard pomerons, 
and employing at larger x the perturbative QCD expressions. The smooth interpolation between 
two regimes for each Q 2 is provided with the help of simple polynomial functions. To obtain 
low-x neutrino-nucleon structure functions F^^ 1 ^ (x, Q 2 ) and singlet part of F^ N,uN (x,Q 2 ) from 
Donnachie-Landshoff function F2 P (x,Q 2 ), we use the Q 2 -dependent ratios R2(Q 2 ) and R^{Q 2 ) 
derived from perturbative QCD calculations. Non-singlet part of F3 at low x, which is very small, 
is taken as power-law extrapolation of perturbative function at larger x. This procedure gives a 
full set of smooth neutrino-nucleon structure functions in the whole range of x and Q 2 at interest. 

Using these structure functions, we have calculated the neutrino-nucleon cross-sections and com- 
pared them with some other cross-sections known in literature. Our cross-sections turn out to be 
the highest among them at the highest energies, which is explained by contribution of the hard 
pomeron. 
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I. INTRODUCTION 



The interest to neutrino-nucleon cross-sections at very high energies, up ~ 10 21 eV, is 
stimulated by High Energy Neutrino Astronomy (HEN A ) (for a review see []], §])■ Ultra High 
Energy (UHE) neutrinos can be of accelerator and non-accelerator origin. In the former case 
UHE protons accelerated in astrophysical sources produce neutrinos in the chain of n- and 
.fT-decays, when UHE protons interact with ambient gas or with low energy photons. Since 
from observations we know that UHE protons exist with energies up to 3 x 10 20 eV [[|, the 
maximum energy of neutrinos is expected up to ~ 10 19 eV. Astrophysical accelerators are 
usually connected with shock waves in SNe, AGNs, GRBs etc, but there could be also some 
other mechanisms of acceleration, such as acceleration in the strong electromagnetic wave 
and in strong electric field due to unipolar inductors (see Ref. [0). 

Non-accelerator sources can provide neutrinos even with higher energies. These sources 
include production by Topological Defects (first suggested in Ref.Q), by decays of super- 
heavy dark matter particles and by annihilation of superheavy particles. Topological De- 
fects in many cases become unstable and decompose to constituent fields, superheavy gauge 
bosons and Higgs particles, which then decay to hadrons and neutrinos. There could be the 
examples when the constituent superheavy fields are produced at annihilation (e.g. annihila- 
tion of monopole-antimonopole connected by string). Annihilation of dark matter particles 
(e.g. neutralino in Earth and Sun) gives another source of high energy neutrino production. 
The maximum energy of neutrinos from above-mentioned sources can reach the GUT scale. 

High energy neutrino radiation from all sources is inevitably accompanied by other radia- 
tions, most notably by high energy gamma rays. Even in cases when high energy photons are 
absorbed in the source, their energy is partly transformed into low energy photon radiation: 
X-rays and thermal radiation. For sources transparent for high energy gamma radiation 
the upper limit on diffuse neutrino flux is imposed by the cascade electromagnetic radiation 
(see Ref. |ll ). Colliding with microwave photons, high energy photons and electrons give rise 
to electromagnetic cascades with most of energy being in the observed 100 MeV - 10 GeV 
energy range. The energy density of this cascade radiation should not exceed, according 
to EGRET observations, u cas ~ (1 — 2) x 10 -6 eV/cm 3 . Introducing the neutrino energy 
density for neutrinos with individual energies higher than E, uj u (> E), it is easy to derive 
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the following chain of inequalities (from left to right): 

A-jr poo 47T 

u cas > uo u {> E) = — El v (E)dE = —EI U {> E). (1) 

C JE C 

An upper bound on integral neutrino flux immediately follows from Eq. ([!]): 

U> E) < (2) 

The latter inequality gives a powerful limit on the possible diffuse neutrino flux. 

UHE neutrinos can be detected at the Earth due to their interactions with nucleons in 
CC and JVC reactions, 

n{n) + iV(e-) -> z=F( e T) + x, (cc) (3) 

i4(*i) + iV(e-) ->!/,(//,) +X, (JVC) (4) 

where I — e,fi, r. These processes also modify the observed neutrino spectrum; neutrinos 
are both absorbed in (|3|) and driven to lower energies in (|4]) on their way from a source to 
a detector ||. 

Cross-sections of z/e-scattering are very small as compared with z/JV-cross-sections. An 
important exclusion is the resonant ^ e e - -scattering pi: 

v e + e~ — > W~ — i> qi + q~j — > hadrons. (5) 

The resonance energy of neutrino is Eq = m^/2m e ~ 6.4 x 10 15 eV; the hadrons are produced 
as a spike with energy Eh = Eq. The number of resonant events in the underground detector 
with the number of electrons N e is given by the simple formula ||: 

v Tes = 2nN e a ef fE I Pe (Eo), (6) 

where <r e // = (37r /y/2)Gp = 3.0 x 10~ 32 cm 2 is the effective cross-section (Gp is the Fermi 
constant) and 2tt is the solid angle open for a deep underground detector (within another 
2tt angle neutrinos are absorbed). 

A detailed discussion of all above mentioned processes can be found in Ref. [[?]]. 

As regards i/JV-cross-sections, especially at extremely high energies, they are unknown 
yet. Really, to calculate the rate of high-energy events in a neutrino detector one actually 
needs the differential cross-sections of u(u)N DIS in the whole range of kinematic variables 
< x < 1 and < Q 2 < oo. Such cross-sections, with QCD-e&ects being taken into 



account (see Ref. ||), may be expressed in terms of Partem Distribution Functions (PDF) 
in proton, in our case quarks, qi(x,Q 2 ), where = u,d,c, s,t,b. However, an influence of 
non-perturbative QCD-eSects on nucleon Structure Functions (SF) at small both Q 2 and 
x cannot be accurately estimated. It makes one to rely just on theoretical models, i.e. on 
various extrapolations. 

In fact, differential neutrino-nucleon cross-section can be parameterized with the help of 
two structure functions, F% N (x, Q 2 ) and F 3 N ' m (x,Q 2 ) (see e.g. Ref.'s g [If], 0). In the 
case of CC scattering @ these cross-sections are 

d 2 af/ N \E V , x, y) _ a S w (1 - y + y ^ N (x, Q 2 ) ±(y- £) x F» N > m (x, Q 2 ) 



dxdy 2 (1 + Swxy) 2 



(7) 



where do = Gp ™ w ; y = Jk, Sw = ^r- and S = 1m^E v , the ± sign corresponds to vjv 
cross-sections, respectively. It is useful to decompose the F 3 N,V into singlet and non-singlet 
parts, 

Ff- N (x, Q 2 ) = F 3 NS (x, Q 2 ) ± F 3 s (x, Q 2 ). (8) 

The NC cross-sections (§) can be presented in a similar way, with mw replaced by mz, Sw 
replaced by Sz = 2mNE v /m 2 Zl a a = G 2 F m 2 z /ii and with structure functions given by 

pvN{NC) g2) = + S 2 + 5 2 + 5l)F „N + {5 2 + ^ _ ^ _ (g) 

if^*, Q 2 ) = (5 2 + 5\ - 51 - 5 2 4 )xF 3 NS . (10) 

The chiral couplings 5i are 

12 11 2 1 

5x = - - -sin 2 w , <5 2 = -- + -sin 2 6W, 5 3 = -- sin 2 6 W , 5 A = - sin 2 9 W ; (11) 

we use sin 2 #w = 0.23117. 



The quark contents of functions FX , F^ 8 and F 3 NS are as follows: 



w 

■ i~ on i~ o /-\t Tiin/ilir\nc £v , ,,,,,, 

2 ) - 1 3 allu - ± 3 

F 2 ^(:r,Q 2 ) = x(q + q), (12) 

^ S (^,Q 2 ) =9-?, (13) 
if (x,Q 2 ) = 2{s-c + b-t), (14) 

with q(x,Q 2 ) = u + d + s + c + b + t and g(x, Q 2 ) — u + d + s + c + b + t. 

Various parameterizations of qi(x, Q 2 ), including recent versions of CTEQ, MRST, GRV, 



may be found in PDFLIB [12|. All these parton distributions were obtained as LO/NLO 



solutions of Dokshitzer-Gribov-Lipatov-Atarelli-Parisi (DGLAP) equations with low Q 2 dis- 
tributions qi(x, Qq) taken from experimental data at Qq ~ 1 GeV 2 . The calculated structure 
functions have been found valid in a wide range of (x, Q 2 ) parameter space: 

1(T 5 <x<l, 1 GeV 2 < Q 2 < 10 8 GeV 2 . 

DGLAP approach is based on perturbative physics. The smallness of the QCD coupling 
constant a s (k±) < 1 implies k\ > Qq. The value of Qq ~ 0.3 GeV determines the smallest 
allowed value of k± and can be viewed as mass of parton in QCD cascade. The minimal 
value of x in perturbative approach is then determined x min ~ Ql/S, where S = 2m^E v 
for z/iV-scattering. 

However, with E v increasing the ever smaller values of x = ^ get important in vN- 
scattering. It should be noted, that HENA actually promises the deepest insight in small- x 
physics. Indeed, the record measurements today by HERA |H| relate to -F 2 6p (a;, Q 2 ) structure 



function with x > 10 -5 , while neutrino-nucleon SF with E v ~ 10 19 eV are sensitive to 
x < 10~ 6 -T- 10~ 8 . As LO/NLO DGLAP dynamics evolves in Q 2 -direction with x being 
fixed, it provides no information about small- £ SF behavior. Moreover, the applicability of 
DGLAP approach to small- a; physics seems to be questionable. It was shown in Ref. WA 



that Mellin transform of DGLAP splitting matrix p(z), p(N) = Jq 1 dzz N p(z), suffers from 
singularities at N = 0. These singularities arise in the perturbative expansion of p q c and 
Pgg in powers of as- 

Nowadays, the only small-x solution obtained in the framework of perturbative QCD 
is the BFKL-pomeron [|16|]. However, the validity of this asymptotic solution is still to 



be checked. Though the importance of this approach is widely recognized, there were 



lately many criticisms of the pure i?Fi^L-pomeron solution (see e.g. Ref. |T^j and references 
therein) . 



A certain modification of i?Fi^L-pomeron description has been proposed in Ref. [pl|. It 
includes a unified BFKLj DGLAP evolution equation with a special 'consistency constraint' 
imposed on BFKL component. Authors applied this approach to calculations of neutrino- 
nucleon cross-sections and checked manifestations of such solution in HENA. 

In this paper we concentrate on the different approach to the small- a; physics, which is 
developing by A. Donnachie and P. V. Landshoff, hereafter DL, with coauthors [111]. It is 



based on the Regge theory inspired description of small-x ep-structure function, -F 2 ep (a;, Q ) 



The authors have actually made the simplest possible assumption, namely, that contributions 
from branch points of the complex /-plane at t — are much weaker than those from poles. 
This hypothesis gives a rather good description of data and may be regarded as a guideline 
in F2 P (x, Q 2 ) small-x extrapolation search. Though, predicting the power-law growth of 
cross-sections at high energies, this approach violates the unitarity. 

According to DL, F2 P (x, Q 2 ) may be written down as a sum of three factorized terms, 

F?(x,Q 2 )=j2f l (Q 2 )x-*% (15) 

i=0 

where eo,ei and e 2 relate to the so-called 'hard pomeron', 'soft pomeron' and g,u,f,a ex- 
changes, respectively. Shapes of fi(Q 2 ) and values of c t are parameters; they are to be chosen 
so that to better fit the data. One of the best and the most convenient for our purposes set 
of/,(Q 2 ),viz. 

^=M^r°H in ( i+ f))' <i6) 

hm - A \w^T^m; (17) 

h(Q 2 ) = A 2 \ (18) 



was proposed in Ref. [14|. With 



e 


= 0.418, 


ei 


= 0.0808, 


C2 


= -0.4525, 


^0 


= 0.0410, 


At 


= 0.387, 


A 2 


= 0.0504, 


«o 


= 7.13, 


a 1 


= 0.684, 


a-i 


= 0.00291, 


X 


= 0.485, 


Ql 


= 10.6, 


Ql 


= 48.0, 



(19) 



it describes 539 data points with x 2 ~ 1-016 per point. 

It is worthy of note, that authors of this approach rather prefer fits with /o oc Q e ° at high 



Q ]l^j. In spite of slightly higher x > such fits are very attractive since their hard pomeron 
term looks similar to the perturbative BFKL solution. However, the direct identification of 
these two pomerons seems to remain dubious. 

In the most recent version of their approach A. Donnachie and P. V. Landshoff 
propose a new parameterization of F2 P (x, Q 2 ) data. They have reduced the number of free 
parameters by including of the real photon cross-section data, cr 7P . Authors argue that this 
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new fit gives a good description of data also for x > 0.001 and for Q 2 < 5000 GeV 2 . With 
no additional parameters, they describe successfully the charm SF, F 2 c (x, Q 2 ), as well. 

However, in the new fit authors use again the power-law dependence of coefficient fo{Q 2 ) 
( pip at high Q 2 . It works good at small Q 2 , which are characteristic for photoproduction 
processes, but it does not extend to Q 2 ~ which are typical for z/iV-scattering at 

S » m^. Such power-law behavior of the hard pomeron term cannot be reconciled with 
the quasi-logarithmic dependence on Q 2 of DGLAP SF at high Q 2 . Therefore we cannot 
directly apply this new successful parameterization to the description of z/iV-cross-sections. 
Maybe it will be possible after a slight modification of the form of hard pomeron coefficient, 
which would coincide at small Q 2 with the original one. But in this paper we are to use the 
previous, logarithmic, parameterization (p^5|)-(p!9|). 

In fact, all descriptions of small-a; SF are basically extrapolations with certain merits and 
demerits. In this paper we construct one more set of structure functions F^, F% N , F% and 
F^ s . We want them 

1. to be defined in the whole range of variables < x < 1 and < Q 2 < oo; 

2. to comprise both perturbative at high-x description, viz. CTEQ5 and non- 
perturbative at low-x pure Regge theory approach of DL; 

3. to be smooth over both variables with limited change of first derivative over logx in 
the interpolation zone. 

Hereafter we call this approach DL+CTEQ5. We hope it is appropriate for the purposes of 
HENA. 

Using DL+CTEQ5 functions, we calculate the CC, NC and total (CC + NC) v{v)N- 
cross-sections and compare them with the results of Ref.'s [0, We also use for comparison 
the cross-sections obtained in the framework of trivial 'logarithmic'' extrapolation, hereafter 
Log+CTEQ5, 

Fr , Lo9+ CTE Q5{x < Q2) = ^CTEQ^^ Q2) ( _X_\^) 5 (2Q) 



d\nx 



%min — 1 X 10 , (21) 



m? n 
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which is analogous to the approach of Ref . . These SF smoothly shoot to the low-x region 
from the CTEQ5 defined high-x one. Starting values of functions and of their logarithmic 
derivatives over x in such extrapolation are taken at the x = x min boundary of CTEQ5. 



II. CONSTRUCTION OF DL+CTEQ5 STRUCTURE FUNCTIONS 

According to (0,|,[l(]) , neutrino-nucleon cross-sections depend just on F% N (x, Q 2 ), 
F^(x,Q 2 ) and F^ s (x,Q 2 ) SF, which quark contents is given by (|I!||14|). However, the 



best small- x data relate to F 2 ep (x, Q 2 )- Being just different combinations of quark density 
distributions, qi(x,Q 2 ), these functions are bound indirectly. So, it is strongly desirable to 
make use of this information. 



A. United F 2 structure function 

Let's first obtain the united, smooth over x and Q 2 , structure function F 2 ep (x,Q 2 ), with 
small-x behavior being in accordance with DL description ( |l5| - ]18f) and with large-x one being 
determined by CTEQ5. These functions are different, of course, so that we are to meet them 
smoothly. 

First, to restrain the Q 2 divergency, we choose the DL parameter set ( |T6| - |T9|) . Really, 
at large Q 2 only the hard pomeron term (|16|) survives in ([15]), so that _F 2 ep (a;, Q 2 ) oc hiQ 2 . 
On the other hand, perturbative dynamics predicts an approximately logarithmic over Q 2 
asymptotic growth of SF as well. In fact, the CTEQ5 Q 2 -dependence of F 2 ep at high Q 2 and 
x pa x min is close to logarithmic, but different; it rather looks like 

F 2 ep (x,Q 2 ) oc (lnQ 2 ) 1+a ^ (22) 

with |«(x)| 1. So, despite both functions fit to the same HERA data at x — x m i n , and 
therefore should coincide in a wide range of Q 2 , they inevitably disperse at large Q 2 . To 
reduce this discrepancy, we solved the equation a(x) =0 numerically; the root is 

x = x ~ 2.527 x 1(T 5 . (23) 

We take the line x(Q 2 ) = xq as one of boundaries of the interpolation zone. It is shown in 
the Fig. |. 



S 
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iog 10 x 

FIG. 1: The borders of interpolation zone corresponding to C = 0.12 (see (|24|)). 

Luckily, xo is rather close to the left CTEQ5 boundary x m i n = 1 x 10 -5 . At x ~ xq 
both DL and CTEQ5 parameterizations still keep valid. It means that at relatively small 
Q 2 < 500 GeV 2 they are to be very close, if not equal. It is important, that at large Q 2 and 
x = xq both descriptions practically do not disperse. 

On having reconciled the Q 2 behaviors, we take concern of smooth SF meeting over x 
at each Q 2 . The x- dependencies of these parameterizations are different. To meet them 
smoothly, we undertake an interpolation between lnF 2 DL and lni 71 2 CT ' E< ^ 5 with the help of 
cubic over In x polynomials; these procedure assures the first derivatives to be continuous in 
the whole interpolation zone. 

So, keeping one border of the interpolation zone at x(Q 2 ) = Xq and varying the shape of 
another border, we call the latter £l,_r(Q 2 ), one gets a set of different interpolations. The 
subscripts L, R mean that at small Q 2 we look for the right border, xr(Q 2 ), while at large 
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iog 10 x 



log i0 Q 2 [GeV 2 ] 



FIG. 2: a) F2 P (x,Q 2 ) as a function of x for different values of log 10 Q 2 = 1,2,3,4,5,6,7,8. Each 
label equals to the value of corresponding log 10 Q 2 - The cubic spline interpolation over x zone lies 
between two dash lines, b) F% p as the function of Q 2 for several values of x denoted in the plot. 
The m\Y line corresponds to Q 2 



m 



Q 2 we look for the left one, xl(Q 2 )- The crossing takes place at Q 2 « 1.78 x 10 6 GeV 2 . 
These borders allow to extend the influence of DL description at small Q 2 to Xq < x < xr, 
on the one hand, and, on the other hand, increase the influence range of the perturbative 
description at high Q 2 to smaller xl < x < x . Such slant border seems us to be reasonable 
from the physical point of view. 

The quality of interpolation may be parameterized by imposing of an additional condition: 



max 



d\nF% N (x,Q 2 ) 



d\nx 



d\nFZ N (x : Q 



<9lnx 



\/{x',x" E[x ,x L:R ]}. (24) 



This actually constrains the maximum change of the tangent in the Xq,Xl_r range. The 
higher is C, the closer borders are allowed. However, at the same time the higher tangent 
change becomes possible in the interpolation zone. An optimum value of this parameter 
seems to be C = 0.12. 

The interpolation zone borders corresponding to this value of C are plotted in the Fig. |l|. 
So, the described procedure results in F 2 ep (x, Q 2 ) with the desired properties. The in- 
dependence of these structure functions is plotted in Fig. 0a for several values of Q 2 . In 
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Fig. 0b the Q 2 -dependence of the same functions is depicted for several values of x. 



B. United neutrino-nucleon structure functions 



Now let us turn to construction of neutrino-nucleon SF F2 N (x,Q 2 ), F^(x,Q 2 ) and 
F^ s (x, Q 2 ). For relation of F^ix, Q 2 ) with F% N (x, Q 2 ) a simple receipt has been proposed 
in Ref. flQfl , hereafter FMR. Under an assumption that z-th flavor quark and anti-quark 



distributions in proton are equal, qi = qi, and that u = d = s = 2c = 2b, the 

^™(x,Q 2 ) = ^Ff(x,Q 2 ) 

rule had been derived there. To get better description, we modify this approach by intro- 
ducing of Q 2 -dependent ratios R2(Q 2 ) and Rz(Q 2 ). These ratio functions may be extracted 
from CTEQ5 at according to the rule 

R2iQ2) = F^ CTE ^(xo,Q 2 ) ' (25) 

MQ2) = F^ TE ^xo,Q 2 ) 5 (26) 

they are plotted in Fig. ||]. These ratios differ from the constant values R2 = and R3 = ^ 
(this value we derived using the assumption of Ref. |10|), denoted in the picture as FMR. 
The difference is especially appreciable at small Q 2 due to the thresholds of heavy quarks 
production. 

Combining ratios ( p5] , p6|) with the constructed -F 2 ep (a;, Q 2 ) and assuming that these rela- 



tions keep valid at lower values of x, we get the small-x 



FZ N (x, Q 2 ) = R 2 (Q 2 ) x F?(x, Q 2 ), (27) 
F 3 s (x,Q 2 ) = i? 3 (Q 2 )xFf(x,Q 2 ). (28) 



At the next step we undertake the analogous to (23) cubic polynomial interpolation. It 



allows to smooth and restrict the x- discontinuity of these parameterizations. 

To describe the negligible small-x non-singlet structure function F^ c (x, Q 2 ), we use 
a trivial extrapolation of the corresponding CTEQ5 function in a way analogous to the 
Log+CTEQ5 (|20|,|2l|) . This completes the construction of united neutrino-nucleon SF. 
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4,5 



FMR 




log 10 Q 2 [GeV 2 ] 

FIG. 3: Ratios R 2 (Q 2 ) = %r and R 3 (Q 2 ) = ^ calculated within CTEQ5 at x = x . The FMR 
lines correspond to analogous coefficients R^ MR = |l and R% MR = A . 

The characteristic features of the derived SF are illustrated in Fig.'s f|a and f|b. 
F% N (x,Q 2 ) are depicted there for several values of x versus Q 2 and for several values of 
Q 2 versus x. 

The behavior of xF%(x,Q 2 ) and of xF^ 1 (x, Q 2 ) as functions of x are demonstrated in 
Fig.'s ^]a and ||b. 

III. COMPARISON OF CROSS-SECTIONS 

Substituting the constructed DL+CTEQ5 SF into Eq.'s (|^J^JTy), we obtain the differen- 
tial i/iV-cross-sections. The following integration over x and y and summation of CC- and 
iVC-inputs gives the total cross-sections as functions of E u . 
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Iog io x log io Q 2 [GeV 2 ] 

FIG. 4: a) F% N {x, Q 2 ) as a function of x for different values of log 10 Q 2 = 1, 2, 3, 4, 5, 6, 7, 8. Each 
label equals to the value of corresponding log 10 Q 2 ■ The cubic spline interpolation over x zone lies 
between two dash lines, b) F% N as the function of Q 2 for several values of x denoted in the plot. 
The mw line corresponds to Q 2 = m 2 ^. 

Denoted as DL+CTEQ5, this sum is shown in Fig. ^. For comparison we also plot here 
the corresponding cross-sections obtained in the framework of 

a) simple Log + CTEQ5 extrapolation (p0|); 

b) CTEQ4 parameterization by Gandhi et al. Ref. M, denoted as GQRS-98 (CTEQ4); 

c) the united BFKL/DGLAP approach by Kwiecinski, Martin and Stasto [113, labelled as 



KMS. 

Due to Regge hard pomeron pole, our approach predicts the more rapid growth of cross- 
sections at high energies. The differences between these calculations become especially clear 
in Fig. 0. We have divided each cross-section by corresponding cross-section of Log+CTEQ5. 
The ratios are plotted in the graph. At E v = 1 x 10 12 GeV the DL+CTEQ5 turns out to 
be twice as high as the logarithmic cross-section. 

This difference is neither unexpected nor dramatic for HENA. One should remember 
that uncertainties in f-fluxes are much higher, while expected low measurement accuracy 
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'°g 10 x lo 9io x 

FIG. 5: a) xF§(x,Q 2 ) as a function of x for different values of log 10 Q 2 = 1, 2, 3, 4, 5, 6, 7, 8. Each 
label equals to the value of corresponding log 10 Q 2 . The cubic spline interpolation over x zone lies 
between two dash lines, b) The same plot as a) but for non-singlet xF^ s (x, Q 2 ). The dashed lines 
correspond to logarithmic extrapolation of corresponding perturbative function as described in the 
text. 

of future detectors and scarce statistics suggest, that such difference may be practically 
insignificant. However, we believe that rapid growth of z/iV-cross-sections may be eventually 
discovered in future giant detectors. This effect may play the essential role for UHE u's 
predicted in the framework of TD models. 

One should also keep in mind, that DL Regge theory approach violates unitarity. It 
implies that the predicted power-law growth of z/iV-cross-sections should be replaced at 
higher energies by, say, a oc In 2 E v one. Though it is yet unknown where and how this 
occurs. 



Conclusions 

In this paper we have derived a new full set of smooth over x and Q 2 ep- and z/iV-structure 
functions, which are defined for arbitrary allowable values of these kinematic variables. 
According to construction, these functions are in agreement with Regge theory inspired 
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FIG. 6: Total (CC+NC) DLGLAP+CTEQ5 z^iV-cross-sections. The analogous cross-sections 



of Log + CTEQ5 extrapolation (20), of CTEQ4 parameterization by Gandhi et al. and of the 
united BFKL/DGLAP approach by Kwiecinski, Martin and Stasto, labelled as KMS are shown 
for comparison. 

hard + soft pomeron small-x parameterization by Donnachie and Landshoff, and coincide 
with perturbative QCD parameterization by CTEQ5 at large x. For the smooth meeting 
of these structure functions over both x and Q 2 , the special interpolation zone boundaries 
have been defined. 

We recalculate the known -F 2 ep (a;, Q 2 ) to (x, Q 2 ) structure functions at small x's with 
the help of introduced ratios R2,z{Q 2 ) ) which are derived from perturbative CTEQ5 descrip- 
tion at x = 2.527 x 1CT 5 . 

Using these new structure functions, we have calculated the i^iV-cross-sections at ex- 
tremely high energies and compared them with those earlier obtained a) within a simple 
logarithmic extrapolation of perturbative structure functions (Log+CTEQ5) and b) in pa- 
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FIG. 7: r{E u ) are ratios of the shown cross-sections and of Log+CTEQ5 cross-section. 



pers |7], 18]. At small and moderate energies these cross-sections are practically indistin- 
guishable. However, at extremely high energies non-perturbative hard pomeron dynamics 
causes a quicker rise of total z/iV-cross-sections with energy. Actually, these growth is the 
highest among all ever predicted in the framework of conventional theories. 

We understand that pure pomeron behavior of SF cannot be a final answer since it 
violates unitarity. Nevertheless, we believe that such approach may be relevant in a wide 
range of energies involved in HEN A. 
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